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ABSTRACT 



Aims. We test by means of numerical simulations and different recipes the efficiency of thermohaline mixing as a process to alter the 
surface abundances in low-mass giant stars. 

Methods. We compute full evolutionary sequences of red giant branch stars close to the luminosity bump by including state of the 
art composition transport prescriptions for the thermohaline mixing regimes. In particular we adopt a self-consistent double-diffusive 
convection theory, that allows to handle the instabilities that arise when thermal and composition gradients compete against each 
other, and a very recent empirically motivated and parameter free asymptotic scaling law for thermohaline composition transport. 
Results. In agreement with previous works, we find that during the red giant stage, a thermohaline instability sets in shortly after 
the hydrogen burning shell (HBS) encounters the chemical discontinuity left behind by the first dredge-up. We also find that the 
thermohaline unstable region, initially appearing at the exterior wing of the HBS, is unable to reach the outer convective envelope, 
with the consequence that no mixing of elements that produces a non-canonical modification of the stellar surface abundances occurs. 
Also in agreement with previous works, we find that by artificially increasing the mixing efficiency of thermohaline regions it is 
possible to connect both unstable regions, thus affecting the photospheric composition. However, we find that in order to reproduce 
the observed abundances of red giant branch stars close to the luminosity bump, thermohaline mixing efficiency has to be artificially 
increased by about 4 orders of magnitude from that predicted by recent 3D numerical simulations of thermohaline convection close 
to astrophysical environments. From this we conclude the chemical abundance anomalies of red giant stars cannot be explained on 
the basis of thermohaline mixing alone. 

Key words, instabilities - stars: abundances - stars: evolution - stars: interiors 



(N 
> 
(N 

m 
oo 

o 
o 



X 



1. Introduction 



After leaving the main sequence, low-mass stars move in the 
HR diagram towards the red giant branch (RGB). During the 
RGB, nuclear reactions take place in a thin shell surrounding 
the helium core and moving outwards in mass. The material 
processed by H-burning is kept hidden inside the core until the 
inner boundary of the convective envelope penetrates deeply 
inwards, reaching the freshly synthetised nucleides. When this 
happens, the material processed by nuclear reactions is dredged 
up to the surface (in the so called first dredge up) modifying 
the photospheric comp osition of red giant stars. Standard 
stellar evolution theory (Ibenl ll967l) predicts no further surface 
abundance variation should take place. However, observational 
evidence strongly suggests the existence of a non-canonical 
mixing processes on the RGB dGilrov 1989I: Gilrov & Brownl 



1991; Lu ckl 11994 ICharbonnell 11994:" 



Charbonnel et al. I1998I 



8; 
0; 



Charbonnel & Do Nascimento 1998; Gratton et al. 2001 
Smith et al.l2002tlShetronel2003HGeisler et al.l2005tlSpite ei"ai] 
2006; lRecio-Blanco & de Lavemvll2007t ISmilianic et al.ll2009l) . 
This extra-mixing seems to be related to the RGB luminosity- 
function bump, i.e., the phase of the evolution when the narrow 
hydrogen burning shell reaches the chemical discontinuity 
caused by the deep penetration of the convective envelope, 
leading to a transitory drop of the luminosity of the star and 
producing a peak in the giant branch luminosity distribution. 



In the last years considerable effort has been devoted to 
identify the non-canonical physical processes that could be re- 
sponsible for modifying the photospheric composition of low- 
mass giant stars at the luminos ity bump stag e . One impor- 
tant clue was first provided by Eggleton et al. (2006) by de- 
tecting the appearence of a mean molecular weight (//) inver- 
sion in a region just above the HBS when the burning shell 
reached the uniform composition layers left behind by the first 
dredge-up phase. Us ing the classic Rayleigh-Taylor criterion, 
Eggl eton et al.l d2006l) found this region to be hydrodynamically 
unstable. The //-inversion detected was identified to come from 
the 3 He( 3 He,2p) 4 He reaction, a process that takes two nuclei and 
transforms them into three, producing a local depression in the 
mean molecular weight per nuclei. This depression is very tiny 
and becomes evident just when it takes place in a background 
of homogeneous chemical composition, like that found by the 
external wing of the HBS at the luminosity bump region. 

ICharbonnel & Zahnl d2007l CZ07) pointed out that, in a star, 
rather than a dynamical instability (Rayleigh-Taylor) it is a 
double-diffusive instability (known in the literature under the 
name of thermohaline instability) what first sets in as the in- 
verse //-gradient builds up. This thermohaline instability takes 
place when the stabilizing agent (heat) diffuses away faster than 
the destabilizing agent (//), leading to a slow mixing process that 
might provide the extra-mixing seeked. 

Since thermohaline instability was identified to take place 
at the luminosity bump several efforts have been conducted 
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in order to understand the actual relevance of this pro- 
cess in modifying the surface ab undance composition of low 
mass giant stars along the RGB (IDenissenkov & Pinsonneault 
2008: ICantiello & Langerl l2010t Iciiarbonnel & Laeardel 12 010: 
Denissen kovll2010T : iDenissenkov & Merrvfi eldll201 ll Stancliffe 
20101) . So, the thermohaline mixing has been studied by means 
of numerical simulations either by considering it as an isolated 
process or combining it with other mechanisms which might 
contribute to the turbulence of the material (e.g., rotation, in- 
ternal magnetic fields). 

Previous investigations are uncertain about the consequences 
and relevance of the setting in of thermohaline instability at the 
RGB: while some authors find the mixing rate generated by this 
process to be enough to reproduce the surface abundances ob- 
served, others encounter this mechanism to be insufficient and 
propose the interaction of more than one process to explain the 
obsevations. Surely a realistic scenario should take into account 
all possible physical processes present and study how they con- 
tribute and interact among themselves. However, there are still 
important doubts in the current treatment of the thermohaline 
mixing as an isolated process which should be addressed before 
we consider actually understood the role of this mechanism. 

One important source of uncertainty comes from the calibra- 
tion of the degree of turbulence generated by each instability. In 
particular, the thermohaline instability gives rise to a slow mix- 
ing of the material which is usually treated as a diffusive process 
characterized by a coefficient that determines the efficiency of 
the mixing. This parameter, the diffusion coefficient, has to be set 
beforehand in order to solve the corresponding diffusion equa- 
tion. To this end, it has u s ually been ad opted a prescription base d 
on the work of lUlrichl (1 19721) and IKippenhahn et all (Il980h . 
where the diffusion coefficient was found to be proportional to 
the square of the (unknown) aspect ratio a (length/width) of 
fluid elements. Unfortunately this means a great uncertainty in 
its value since the linear theory does not give a reliable estimate 
of the maximum length of salt fingers relative to their diameter. 

Being the subject still a matter of debat e, laboratory exper- 
iments simulating oceanic conditions (e.g., Krishnamurti 2003) 
suggest a geometry of slender fingers for the convective elements 
and thus some authors adopt high values of a (> 5) in order to 
reproduce the surface abundances of low-mass stars after the lu- 
minosity bump. However, physical conditions inside a star are 
very different from those in the laboratory and it is not clear if 
elongated structures can be stable, especially when shear and 
horizontal turbulenc e is present. In view of these concerns, other 
authors (e.g.. IKippenhahn et al.ll 19801: ICantiello & Langerl2010l) 
adopt "blobs" (ax 1) as the preferred fluid element morphol- 
ogy. This freedom in the election of the aspect ratio has an evi- 
dent impact on the diffusio n coefficient which has be en reported 
to affect the results (CZ07. ICantiello & Langerll2010l) . 

Very recently, IDenissenkov! (1201 Ol) and iTraxler et aD (1201 ll) 
have presented the first numerical simulations of thermohaline 
(fingering) convection close to the astrophysical regime — i.e. 
Prandt l number Pr ~ 10~ 6 and inverse Lewis number t ~ 10~ 6 . 
In fact, Traxle r et al ] (l20Tl .bv means of high performance three 
dimensional simulations, were able to derive asymptotic scaling 
laws for thermohaline composition transport. These asymptotic 
scaling laws are, then, the first empirically motivated and pa- 
rameter free available recipe for the treat ment of therm o haline 
mixi ng in an ast r ophys ical regime. Both IDenis senkov (2010) 
and ITraxler et al.l (1201 ll) have suggested that their results imply 
that the thermohaline mixing is not efficient enough to account 
for the changes in the surface abundances of red giants close to 
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Fig. 1. Regions of different stability regimes for diffusion rates 
A — 0, D = 0.001 and F — 0. Note the convective region with 
V - V a( j < for which the standard mixing length approach 
(Vp = 0) can not provide mixing velocities and that have been 
sometimes misidentified with the thermohaline regime. 



the luminosity bump, but no stellar evolution computations have 
been performed. 

In the present work we test these suggestions by means of 
full evolutionary simulations of the development of thermoha- 
line convection in RGB stars. Specifically we study the rel- 
evance of thermohaline mixing in RGB stars when more so- 
phisticated and physical ly sounding prescriptions, than that of 
Kippenhahnet al. (1980), are adopted in a stellar evolutionary 
code. In par t icular we adopt the very recent prescription of 
ITraxler et al.l (1201 1[) and th e double diffusive mixing length the- 
ory of Grossman & Taaml (Il996l) . While the former is based on 
realistic 3D numerical experiments and is essentially parameter 
free, the latter successfully reproduces most previously known 
results about convection in astrophysics and, when composition 
gradients are considered, it establishes its own stability condi- 
tions, thus providing a new perspective to study the thermohaline 
instability problem. 



2. GNA convection theory 

As an effort to p rovide a better nonlocal theory of convection, 
iGrossman et al.1 (1 19931) developed a flexible and powerful for- 
malism, which was designed to make unbiased, self-consistent 
predictions about complex phenomena associated to the trans- 
port of energy in stars. Here we use this fo rmalism and follow 
the prescription of Grossman & Taaml (fl 996) to get the local the- 
ory of convection in a composition-stratified fluid. 

Basically the theory allows to find the mixing rate of the fluid 
in the convective, thermohaline and semiconvective regimes by 
solving two equations simultaneously: the first of these equation 
corresponds to the turbulent velocity cr, 



cr- 



(A+D + 2Bo-)— (V - V ad ) - (A + F + 2Bcr)— V 
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Fig. 2. Profiles of the abundances of H, 3 He, 12 C, 13 C, 14 N, 16 
and of the mean molecular weight gradient as a function of 
mass coordinate. The full line stands for V^ > and the plus 
signs otherwise. Left and right panels correspond to the situation 
before and after the bump, respectively, for a Q.9Mq model with 
1 x 1CT 3 . The abscissa ranges from the bottom of the hydrogen 
burning shell to the base of the convective envelope. 



-(A + D + 2Bcr){A + F + 2Bcr)(D + F + 2Bcr) 



(F + Bcr)§^(V - V ad ) - (D + B&)§-V„ 



-{A + Bcr)(D + Bo-)(F + Bcr) 



= 0, 



(1) 



where V = dlnT/d\nP, V ad = (3 In T/d In P) ad is the adiabatic 
gradient, V^, = d Inp/d In P is the molecular weight gradient, g is 
the local acceleration due to gravity, a = -{d In p/d In T)p tfl is the 
coefficient of thermal expansion, <p = (d\np/dlnp)pj, H v is the 
pressure scaleheight, D, F and A are the diffusion rates of heat, 
composition and momentum, respectively, and B = 2/1, with 
/ the u nique mixing length considered by iGrossman & Taarnl 
( 1996). The other equation involves the flux conservation 



v Rad - v ad = (V - v ad ) + h v (^pi we, 



(2) 



where VR ad is the temperature gradient that would be necessary 
to transport the whole flux by radiation, p is the density, Cp is 
the specific heat, K is the radiative diffusive conductivity, and 
w8 is the correlation between turbulent velocity and turbulent 
temperature excess given by 



( * 

— - \ D+F+lBo- 
W0 = 7 



. A - F _2Bcr- ¥ ^)(A + Bo-)-Y Ta i 



D+F+2Bcr 



-A-F- 2Bct - 



+ Y 



ag 



(3) 



where, for simplicity, we have set that X = ga(V - V ad )/// p and 
Y = gipV^/Hj,. Our Eq. (0 differs from Grossman & Taam's 
Eq. (17) because we have fixed some sign errors present in the 



original expression. It is worth mentioning that for realistic cases 
D » A, D » F and D goes to zero. 

Equations ([]]) and © have to be solved simultaneously for cr 
and (V - V^). The leading factor of Eq. ([TJ shows that cr = is 
always a solution, but it corresponds to a stable equilibrium only 
if no other real and non-negative root exists. In general, the fluid 
will seek out the most turbulent equilibrium state, thus, if more 
than one root is positive the system will evolve to the largest 
root, being the cr = solution unstable. 

In the present work the mixing of nuclear species of mass 
fraction X, is performed by solving the diffusion equation 



dXj 
dt 



(dX, 


| ! d 


\ dt 


L + Mr 



(4nr 2 p) 2 D c ^- 
v H ' dM r 



(4) 



with the diffusion coefficient D c define d in terms of the tur bulent 
velocity cr and the mixing length I by dWeaver et al.ll 1978b 



1 



-crl. 



(5) 



Appendix A contains a few additional details about the proce- 
dure followed by us in order to solve GNA's equations. 

3. An empirical scaling law for compositional 
transport by fingering convection 

Although double-diffusive processes have been studied 
by several authors by me a ns of hydro d ynamics codes 
(see, e.g.. iMerrvfieldl Il995t iBiel ol 1200 It iBascoull 12007b 
IZaussinger & Spruitl 1201 Ol) . It was 



Traxler et all (1201 lb who 



performed the first three-dimensional simulations to address the 
question of double-diffusive transport by fingering convection 
in astrophysics. It is important to note that Trax ler et al.l (1201 ll) 
conducted their simulations at Pr ~ O(10~ 2 ) while the true 
astrophysical regime occurs at Pr ~ O(10~ 6 ). Therefore, their 
empirical scaling law relies on the validity of the asymptotic be- 
havior suggested by their results. They model a finger-unstable 
region using a local Cartesian frame (x,y,z) oriented so that its 
vertical axis z has a direction opposite to that of the gravitational 
acceleration. Also the Boussinesq approximation is used and, 
consequently, it is assumed that small density, temperature 
and compositional perturbations (p, T, jX) are related by the 
following linearized equation 



— = aT + Bp., 
Po 



(6) 



where po is a reference density, a = -p Q l dp/dT, and B = 
p^dp/dfx. Expressing the velocity, temperature and composi- 
tional fields as a background component plus a perturbation, it is 
obtained 



u(x,y,z,t) = u(x,y,z,t), 
T(x,y,z,t) = T (z) + T(x,y,z,t), 
H(x, y, z, f) = fioiz) + p(x, y, z, t), 



(7) 
(8) 
(9) 



with Tq(z) - zdT/dz and /jq = zd/j/dz. By scaling the time 
(f), the temperature and the c omposition adequat ely by means of 
the expected finger scale (see lTraxler et al.l201 It for details), the 
final set of equations to solve turns out to be 



1 da . 
— — + u • Vu 

Prl dt 



-Vp + (f -p)e z + V 2 u 



(10) 
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V • u = 0, 



— + w + u ■ Vf = V 2 f, 
dt 

du w _ —2. 

+ — + u • Vu = tVV 
dr /? 



(11) 
(12) 

(13) 



where w is the z component of u, Pr is the Prandtl number, p 
is the non-dimensional pressure perturbation from hydrostatic 
equilibrium, Rq = (V - V ac j)/V jU and r = k^/kj, with the com- 
positional diffusivity (see below) and kj the thermal diffusivity 
given by 



4acT i 
3kCpp 2 



(14) 



In this last equation a stands for the radiation density constant, c 
is t he speed of light and k the Rosseland mean opacity. 

iTraxler et alJ (1201 1 ) solved Eqs. (I10til3b in a triply periodic 
box of size (L x ,L y ,L z ) and carried out simulations for moder- 
ately low values of the Prandtl number and diffusivity ratio of 
order 0{\Q 2 ). 

As a result of numerical experiments, it turns out that the tur- 
bulent compositional transport by fingering convection follows a 
simple law for the diffusion coefficient, namely 



D„ 



101 V^ve" 3 - 6r (l 



r) 1 



(15) 



where r - (Rq - 1)/(t 1 - 1) and v is the total visc osity given by 
the su m of the molecular and radiative viscosities (iDenissenkovl 
12011 



V = V mo i + v rad 
with 

4aT 4 



Vrad = 



and 



15 CKp 2 



(16) 



(17) 



J5/2 



vmoi = k. = 1.84 x 10 _17 (1 + IX)—, [cmV 1 ] (18) 

P 

where X is the hydrogen mass f raction. Based on the asymptotic 
behavior shown by their results. ITraxler et aLl(l201 lb suggest the 
possibility of applying Eq. (TT3T > to the more extreme astrophysi- 
cal regime, provided Pr is of the order of r. 



4. Numerical simulations 

In order to study the effects of thermohaline instability in 
low-mass giant stars we performed si mulations using a o ne- 
dimensional evolution code (LPCODE, lAlthaus et al.ll2005l) in- 
corporating GNA's convection theory to compute the mixing 
rates of the different stability regimes defined by this formalism. 
We adopted in our numerical experiments the following choice 
for the parameters: A = 0, F = 0, D = 3K/(pC P l 2 ), a = 1, 
<p — 1 and I = 1.35 (approxim ately equivalent to a mixing l ength 
parameter of 1.61 in the usua 1 Kip penhahn & Weigertt 1990 pre- 
scription), and impleme nted the same nuclea r reaction network 
as described at length bv lAlthaus et all (120051) . 

We computed stellar models of 0.9 Mq, 1.3 Mq and 1.6 
Mq, each one with three different initial metallicities, namely 
Z=3.17 x 10~ 4 , 1 x 10~ 3 and 6.32x 10~ 3 , and let them evolve from 
the main sequence until after the luminosity bump. For each 



model, we paid special attention to the detailed stellar structure 
of the region between the HBS and the base of the convective en- 
velope, contrasting the situation given before the star enters the 
luminosity bump region and after that stage of its evolution. The 
results obtained in all cases were qualitatively very similar, thus 
we show here just one case, namely the 0.9 Mq, Z=l x 10~ 3 
model, which is representative of what happens to the others. 
Figure|2]shows the abundance profiles of some elements and the 
mean molecular weight gradient for the 0.9 Mq, Z=l x 10~ 3 
model, before and after the luminosity bump. Before the lumi- 
nosity bump (left panel in Fig. the mean molecular weight 
gradient shows two peaks, corresponding to the hydrogen 
burning shell (left peak) approaching the molecular weight dis- 
continuity (right peak) left behind by the first dredge-up. When 
the hydrogen burning shell reaches the discontinuity, the reac- 
tion 3 He( 3 He,2p) 4 He produces a molecular weight inversion in 
the external tail of the HBS, destabilizing the region and produc- 
ing thermohaline convection. The destabilized zone (shown by 
plus signs in the right panel of Fig. 13 never reaches the convec- 
tive envelope in our simulations, being both regions separated by 
a radiative zone that prevents any change in the surface compo- 
sition of the star. This result clearly differs from the main result 
presen ted by CZ07 when the slender finger geometry of Ulrich 
dl97l was adopted. This should not come as a surprise as our 
assumption of a unique mixing length in GNA theory is far from 
a slender finger geometry. In fact, our results are consistent with 
those of CZ07 when blobs, rather than slender fingers, are as- 
sumed. As shown by CZ07, different blob/finger geometries can 
affect the diffusion coefficient by more than two orders of mag- 
nitude. In this connection, we performed additional simulations 
artificially increasing GNA's diffusion coefficient of thermoha- 
line unstable layers by a factor of 10 3 in order to test the even- 
tual relation between the more rapid mixing rate and the surface 
abundance variations. Figure [3] shows the abundance profile of 
the same elements included in Fig. [2] as well as the run of the 
molecular weight gradient in the region comprised by the HBS 
and the base of the convective envelope, for this new experi- 
ment. Now the thermohaline zone expands outwards (in mass) 
occupying all the former radiative region that separated it from 
the convective envelope. The contact between both convective 
regions allows for the non-canonical extra mixing to take place, 
thus modifying the photospheric chemical composition after the 
luminosity bump. 

Given the contrast between the diffusion coefficients com- 
puted as given by the GNA and those reported by CZ07, we 
decided to perform further simulations combining the GNA the- 
ory with other prescriptions used to estimate the diffusion coef- 
ficient. Only very few different prescriptions for the computa- 
tion of the diffusion co efficien t in the rmohaline unstable regions 
exist in the literature. lUlrichl (11972b was the first to derive an 
expression f or the turbulent diffusivity produced by that instabil- 
ity, whereas iKippenhahn et alJ dl980h extended previous works 
to the case of a non-perfect gas. The linear theory used by these 
early works yield a solution for the diffusion coefficient that is 
proportional to the square of the unknown aspect ratio (length to 
diameter) of the fluid elements, an issue that is still a matter of 
debate. Indeed, the diffusion coefficients may differ by about two 
orders of magnitude depending on the form factor adopted by 
different authors. Thus, the implementation of the linear theory 
has the drawback of containin g a high intrinsic u ncertainty. On 
the other hand, more recently. ITraxler et al.l (1201 lb successfully 
derived empirically determinated transport laws for thermoha- 
line unstable regions by means of three-dimensional simulations 
performed at parameter values approaching those relevant for 
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Fig. 3. Same as Fig.[2]for an artificially increased diffusion coef- 
ficient (see text). 

astrophysics. This represents an alternative and more physically 
sound approach which helps us to avoid the problems of the clas- 
sical linear theory and supplies an independent way to address 
the question of the actual role of mixing in thermohaline unsta- 
ble regions. 

Since GNA convection theory may be implemented to deter- 
mine the regime of energy transport of any layer with the ad- 
vantage of leaving the computation of the diffusion coefficient 
as an independent task which might subscribe to different pre- 
scriptions, we decided to study the system's response combining 
GNA formalism with two independent recipes. On the one hand, 
we computed the thermoh aline diffusion coefficient by means of 
the expression obtained by Kippenhahn et al. ( 1980) 



D K = a t h 



3K 



2pC P (V-V ad )' 



(19) 



where a± is a efficiency parameter which depends on the ge- 
ometry of the fluid elements, p is the density, K = 4acT 3 /(3/cp) 
the thermal conductivity, and Cp = (dq/dT)p the specific heat 
capacity. We s et q t h = 2, which rough ly corresponds to the 
prescription of iKippenhahn et al.l d!980l) . On th e other hand, 
we co mputed diffusion coefficients adopting the iTraxler et al.l 
(1201 ll) empirical law given by Eq. (fl3T l. 

Thus, we performed a few additional simulations for the 0.9 
Mq, Z=l x 10~ 3 , and 1 Mq, Z=0.02, sequences in order to inves- 
tigate the response of the system when we solely vary the recipe 
to compute diffusion coefficients. Fig. [4] shows the evolution of 
th e thermohaline reg i on alo ng the RGB when the prescription 
of Kippenhahn et alJ dl980h is adopted. Note that the convec- 
tive envelope never enters into contact with the thermohaline 
region. Consequently, for this model and mixing treatment, the 
photospheric abundances of the star remain constant through- 
out this phase. A similar behavior is shown by Fig. [5] corre- 
spond ing to the implementation of the recipe of Traxle r et alJ 
(1201 ll) . In this case, the thermohaline zone is much more nar- 
row than before. We will see in the next section that this fact is 
in close relation with the magnitude of the diffusion coefficients 
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Fig. 4. Evolution of the region between the HBS and the convec- 
tive envelope, when the prescription of Kip penhahn et al.l |l980) 
is adopted to compute the diffusion coefficient in the thermo- 
haline zone. Time interval spans from the instant when the stars 
luminosity reaches L m 96 (i.e., before the luminosity bump) un- 
til L m 1826, close to the top of the RGB. The figure corresponds 
to the 0.9 Mq, Z=l x 10~ 3 , model. 



computed using different prescriptions. Finally, it is worth not- 
ing that in our 1 Mq , Z=0.02, sequence, we did not find any 
contact between the bottom of the convective envelope and the 
thermohaline regio n. This result is at varian ce with the simula- 
tions presented bv lCantiello & Langerl(l2010l) which showed that 
this co ntact occurred in 1 M^ m ass stars even for the prescrip- 
tion of Kippenhahn et alJ dl980h with a± — 2. We suspect that 
this different behavior may be due to the different microphysics 
assumed in both stellar codes. 



5. Summary and discussion 

We have studied the impact of thermohaline mixing in red giants 
close to the luminosity bump in the light of two non-standard 
and physically soun ding mixing prescriptions: the GNA and the 
Trax ler et al.l (1201 lb prescription. To the best of our knowledge, 
this is the first time that the empir i cally b ased thermohaline mix- 
ing prescription of Tra xler et al.l (1201 lb is tested in the context 
of detailed evolutionary simulati ons. In the case of the double 
diffusive mixing length theory of [Grossman et al] (1 19931) it has 
allowed us to include thermohaline mixing in a consistent way 
with the other unstable regimes that are possible when + 0, 
selfconsistently solving the temperature gradients and turbulent 
mixing rates. 

For the sake of completeness let us mention that in the case 
of GNA theory, we fi nd the thermohaline mixin g efficiency to be 
very similar to that of Kippenhahn et al .(1980). In fact our com- 
putations show tha t at al most all layers the value predicted by 
Grossman & Taam1(ll996l) is £>gna ~ ^Kip/6, a difference that is 
just a consequence of different choices in the adimensional coef- 
ficients of both prescriptions. The similarities between these two 
prescriptions should not come as a surprise as the GNA theory 
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Fig. 5. Same as Fig.[4]but when the prescription of iTraxler et al.l 

Ml is adopted to compute the diffusion coefficient in the ther- Fi §- 6 - Relatlon between diffusion coefficients computed using 

mohaline zone. The figure corresponds to the 0.9 M Q , Z=l x the_prescnptions of Kippenhahn et al. ( 1980), £> Kip ,| Traxler et al. | 

10 -3 model si (1201 II) . Dt„. and lGrossman et al.l (119931) . ZV.ma . 



is a sophisticated version of the mixing length theory but still 
relies on a very similar p icture than the standard MLT, in which 
Kippenhahn et al. ( 1980) prescription is based on. We consider 
our results as an actual validation o f the GNA theory for the 
cases where Kippenhah n et al.l dl980l) prescription i s applic able. 

Both iTraxleret al.l (1201 ll) and iGrossman etaD (119931) pre- 
scriptions have identified thermohaline mixing to develop in 
RGB stars close to the luminosity bump, in agreement with all 
previous work that have adopted more simplified a pproaches 
dCharbonnel & Zahnll2007t ICantiello & Lang er 2010). 

However , in a greement with iDenissenkovl (1201 Ol) and 
ITraxler et al.l (1201 lb suggestions, our full evolutionary calcula- 
tions confirm that thermohaline mixing is not efficient enough 
for fingering convection to reach the bottom of the convec- 
tive envelope of red giants. Thus, no changes in the surface 
chemical abundanc es of red giants are obt ained when either 
ITraxler et al.1 (1201 ll) or|Grossman et al.l (Il993l) prescriptions are 
adopted. Interestingly enough, as the value of (V - V a d)/V /J in 
the thermohaline zone is (V - V a d)/V /J ~ 10 3 .. . ~ 10 4 it falls in a 
regime in which the standard prescription of Kippen hahn et alJ 
J 19801) strongly over estimates the therm ohaline mixing effi- 
ciency (see Fig. 3 of Tra xler et al.l 1201 ll) . As can be seen in 
Fig. (O the standard prescri ption is ~ 100 to 1 000 times more 
efficient than t he empirical ITraxler et al.l (1201 ll) law. However, 
we know from lCantiello & Langer! (1201 Ol) that the standard pre- 
scription is still not enough to account for the surface abun- 
dances of RGBs. Thus, in order to allow contact between the 
thermohaline region and the convective envelope, the diffusion 
coefficient should be about 4 orders of magnitud e higher than 
predi c ted by realistic ther mohaline transport laws (IDenissenkovl 
l2010t ITraxler etldl 1201 ll) . Since hydrodynamics codes have 
shown to be consistent, yielding similar results between dif- 
ferent implementations, they should be trusted in the physi- 
cal regime studied (Pr > 10 2 ) which, due to computational 
limitations, is not the actual astrophysical regime (Pr~ 10 6 ). 



While the prescriptions used here still relies on an asymp- 
totic scaling, it seems difficult that the diffusion coefficients 
could be off by this much. Thus, we can conclude that thermo- 
haline mixing alone is very unlikely to be the explanation for 
the chemical abundance anomalies of red giants. 
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Appendix A: Solving GNA's equations 

GNA theory of convection provides us with a set of equations 
that have to be solved in order to find the values of the tem- 
perature gradient V and the turbulent velocity cr of the system 
in regions of different energy transport regime. In practice, this 
means to solve equations ([]]) and (f2| simultaneously. Being it im- 
possible to express any of these variables in terms of the other we 
followed an iterative procedure, adopting the Newton-Raphson 
method to this end. In order to avoid eventual numerical insta- 
bilities associated to the divergence of Eq. Q when the denom- 
inator becomes small, we elementary transformed equation (f2]l 
by multiplying it by that denominator, and rearranging the flux 
conservation equation we obtain 



a\ and «2 are the coefficients of the quadratic equation in X (i.e., 
V), which depend explicitly on the turbulent velocity cr, the com- 
position gradient V^, and the total radiation gradient VR a d. 

Thus, given a set of diffusion rates of heat (D), composition 
(F) and momentum (A), and once and Y are known, we 
first determine if the total radiation might be transported in a 
not convective way. If radiative transport is insufficient, convec- 
tion has to carry some fraction of the energy flux, and thus we 
start the iterative procedure above mentioned. We adopt an ini- 
tial (guess) value for X and solve Eq. (Q]i for cr. As stated before, 
the system will seek out the most turbulent equilibrium state, 
thus we solve both cubic equations and pick up the largest pos- 
itive root. The adopted values for X and cr are then introduced 
in Eq. (1A.U and Newton-Raphson method is used in order to 
find the correction to be applied to X. Iterating this process it is 
possible to obtain the values of X and cr that satisfy equations 
(Q} and (lA.ll i. Numerical experiments have shown that X - X^d 
is a good starting value for the Newton-Raphson process, while 
other choices resulted in false roots found by the algorithm. 

Finally, it is worth mentioning that factors in brackets in Eq. 
(Q]l are cubic in cr and, since both cubics are different, the condi- 
tions that separate the real roots region from the one real plus two 
complex conjugate roots region are also different. Despite this 
difference, for the stellar astrophysics case we have that D » A 
and D » F, and both conditions tend to the same curve, thus 
being unnecessary in practice to compute both limiting curves. 



X 1 + fll (cr, Y, X R . dd ) X + a 2 {(T, Y, X Rad ) = 0, 



(A.l) 



where we adopted the following nomenclature 



X = ga(V - V ad )/# p , 

XRad = gC^VRad - V ad )/i/ p , 

Y = g^/Hp. 



(A.2) 
(A.3) 
(A.4) 



